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Abstract 

Low-energy orbital magnetic dipole excitations, known as scissors mode (SM), are studied in 
alkali metal clusters. Subsequent dynamic and static effects are explored. The treatment is based 
on a self-consistent microscopic approach using the jellium approximation for the ionic background 
and the Kohn-Sham mean field for the electrons. The microscopic origin of SM and its main features 
(structure of the mode in light and medium clusters, separation into low- and high-energy plasmons, 
coupling high-energy Ml scissors and E2 quadrupole plasmons, contributions of shape isomers, etc) 
are discussed. The scissors Ml strength acquires large values with increasing cluster size. The mode 
is responsible for the van Vleck paramagnetism of spin-saturated clusters. Quantum shell effects 
induce a fragile interplay between Langevin diamagnetism and van Vleck paramagnetism and lead 
to a remarkable dia-para anisotropy in magnetic susceptibility of particular light clusters. Finally, 
several routes for observing the SM experimentally are discussed. 
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I. INTRODUCTION 



This paper deals with orbital magnetism in metal clusters. Because of a possibly large 
number of atoms, the valence electrons may accede single-particle orbitals with very high 
angular momenta. The occupation of these orbitals has a large impact in cluster static 
magnetism |2], [3] or in collective magnetic modes of orbital nature [|, ^ ||. Two 
remarkable examples are the scissors |1|, f|, |5| and twist |6[] modes. The scissors mode is 
strictly correlated with cluster deformation. It can be viewed as a small- amplitude rotational 
oscillation of a spheroid of valence electrons against a spheroid of the ionic background (hence 
the name SM). 

The SM is a general dynamical phenomenon already found or predicted in different quan- 
tum finite systems. It was first proposed [0] and observed |§ in atomic nuclei where it still 
remains a hot topic for both experimental and theoretical studies (for a review see ||). 
It was later predicted in a variety of different systems, like metal clusters |], f|, quantum 



dots |L(J and ultra-cold superfluid gas of fermionic atoms |TlJ. More remarkably, it was 
predicted |L2] and observed |U| in a Bose-Einstein condensate. All these different systems 
have two features in common: The broken spherical symmetry and the two-component na- 
ture (neutrons and protons in nuclei, valence electrons and ions in atomic clusters, electrons 
and surrounding media in quantum dots, atoms and the trap in dilute Fermi gas and Bose 
condensate) . 

The SM strongly affects orbital magnetism in alkali metal clusters. It can be described 
in terms of Larmor diamagnetism and temperature-independent van Vleck paramagnetism 
|l|, 0. They are both weak and, therefore, need to be studied in systems, like alkali metal clus- 
ters, where strong forms of magnetism, like ferromagnetism, are absent. Being a low-energy 
mode, the SM determines the van Vleck paramagnetism and causes a strong anisotropy 
in the magnetic susceptibility |], 14 1. Moreover, some particular light clusters can exhibit 
dia-para anisotropy, being diamagnetic along the z symmetry axis and paramagnetic in x-, 



y-directions [[14 



The SM has already been studied in schematic |jj and microscopic ||, [14], [L5|] approaches. 
The microscopic calculations, though accounting for quantum shells effects, were not fully 
self-consistent. A deformed Woods-Saxon fll6| , rather than a self-consistent Kohn-Sham, 
one-body potential was adopted. The quadrupole deformation was deduced from other 
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models or experimental estimates. Certainly, we need to perform fully self-consistent calcu- 
lations based on density-functional theory in order to settle the subtle issues, like the fragile 
dia-para anisotropy, the role of the ionic structure, triaxiality, shape isomers, etc.. 

Self-consistent calculations, accounting for the ionic structure, were performed for spin- 
saturated ground and spin-polarized isomeric states of light sodium clusters Nai2 and Nai6 
H . It was found that the scissors response remains determined basically by the global defor- 
mation, in spite of the fact that triaxiality and ionic structure induce a strong fragmentation 
in the strength. It was also shown that the detailed ionic structure destroys locally spher- 
ical symmetry thereby causing a finite, though very weak, Ml response (transverse optical 
mode) even in clusters with zero global deformation. 

Calculations of the same level of completeness become prohibitive as one moves to heav- 
ier clusters. For this reason, we are forced to use in the present paper the Kohn-Sham 



approach with a soft jellium model for the ionic background ^L7j. This simplifies greatly 
the calculations and allows to proceed to heavier clusters. At the same time, the jellium 
approach is accurate enough for the principle problems considered here. The treatment of 
the electrons is fully self-consistent. We adopt a deformed Kohn-Sham mean field using the 
energy functional of ||18|| . The cluster shape (in terms of axial quadrupole and hexadecapole 
deformations) is determined by varying the jellium deformation and minimizing the total 
energy of the system. The optical response in the linear regime is calculated within the 
separable random-phase-approximation (SRPA) method self-consistently derived 

from the Kohn-Sham functional. The SRPA has been already successfully employed for the 
description of the dipole plasmon in spherical and deformed p0|, f22| alkali metal clusters. 

In the present paper, we will consider both the optical Ml response and static magnetic 
orbital effects. In Section 2, macroscopic and microscopic treatments of the SM are briefly 
outlined. In Section 3 the calculation scheme is presented. The optical response is discussed 
in Sections 4 and 5. It will be shown that, in analogy with the electric dipole plasmon in de- 



formed clusters |2(| p2]j , the Ml response in light clusters has a distinctive profile determined 
by the deformation. Instead, due to a strong Landau damping and contributions of shape 
isomers, the response in medium clusters becomes vague and the SM can be viewed as a 
statistical mix of contributions from different cluster shapes. Besides, we discuss structure of 
the residual interaction and explain a small collective shift in SM excitations. The coupling 
between the high-energy SM branch and the electric quadrupole plasmon is demonstrated. 
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In Section 6, the dia-para anisotropy is discussed. In Section 7, we estimate perspectives 
to observe the SM in photo-absorption, Raman scattering and inelastic electron scattering. 
The conclusions are given in Section 8. 

II. THE SCISSORS MODE: BRIEF OUTLINE 

The macroscopic and microscopic treatment of the SM in clusters is discussed in detail 
in Refs. Jl|, f|, ||. Thus, we give here only a brief outline for a better understanding of the 
results presented in the next sections. 

In the geometrical model of J7|, the SM arises from a rotational oscillations of valence 
electrons versus the ions, both assumed to form distinct spheroids (see left part of Fig. 1). 
Following the alternative view of |1| (right part of Fig. 1, the displacement field of the mode 
is a sum of the rigid rotation and a quadrupole term (the latter provides vanishing velocity 
of electrons at the surface): 

u(r) =Q x r + 6 2 (l + 6 2 /3)- 1 V(yz) (1) 

where 5 2 is the quadrupole deformation parameter (to be defined in the next section). In 
axially symmetric systems, the SM is generated by the orbital momentum fields L x and 
L y perpendicular to the symmetry axis z, and it is characterized by the quantum numbers 
| A 71 " = 1 + > where A is the eigenvalue of L z and 7r is the space parity. Energy and magnetic 
strength of the mode can be estimated macroscopically [p]. ||: 

u = ™N?'% eV, (2) 

B(M1) = 4(1+ | L x | 0)V' 

= \n £ <r 2 > Uf4 (3) 
- N^h Hi 

where N e is the number of valence electrons and r s the Wigner-Seitz radius. We use here 
natural units m e = % = c = 1. The value B(M1) stands for summed strength of the 
degenerated x- and y-branches. The z-branch vanishes for symmetry reasons. It is worth 
noting that B(M1) does not depend on r s and so is the same for different metals. 
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The microscopic treatment of the SM yields the shell structure of an axially deformed 
mean field. One can characterize the emerging single-electron states in terms of the quantum 
numbers of the axially deformed harmonic oscillator (Clemenger-Nilsson basis). These are 
the triplets v = [Afn z A] where n z labels the number of nodes in z-direction (=symmetry 
axis) and Af is the principle shell number Af = n z + 2n r + A (from which one can derive 
the number n r of radial nodes). The angular momenta orthogonal to the symmetry axis, L x 
and L y , promote low-energy AAf = transitions inside the valence shell and high-energy 
AAf = 2 transitions across two shells. Moreover, one may expand wave functions in terms 
of the spherical basis (nLA) 

^ =[A/ -„ Z A] = < L RnL(r)Y LA (n). (4) 

nL 

This allows to evaluate the single-particle orbital Ml transition amplitude between hole 
[y = h) and particle {y = p) states 

(%\L x \V h ) oc 5 w 5 AptAKii (5) 
£ <L*nL y/L(L + l)-A h (A h ±l). 

nL 

Eq. (^) shows that the scissors mode is generated by A p = ± 1 transitions between 
the components of one and the same spherical (nL)-level. In spherical systems, (nLA)- 
states belonging to the level {nL) are degenerate while in deformed systems they exhibit the 
deformation splitting and so may be connected by Ml transitions with non-zero excitation 
energies. This is the origin of the scissors mode. The energy scale of the scissors mode 
is determined by the deformation energy splitting and so is rather small. This explains 
the predominantly low-energy (AA/"=0) character of the scissors mode. Just the low-energy 



branch carries most of the scissors B(M1) strength (see also discussion in Section |TV|). The 
high-energy (AAf =2) branch of the mode is much weaker since the particle states involved 
into (AAf =2) transitions include only small (nL)-components from the valence shell. 



III. CALCULATION SCHEME 



Our approach |2(J employs the Kohn-Sham equations for the electronic mean field using 



actually the energy-density functional of [18|. The positive ionic background is modeled by 
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TABLE I: Ground state deformation parameters 5 2 and £4 and moments j3 2 an d /?4- For Na + 5 and 
Na + 19 the deformation parameters for the isomeric states together with their energy deficits AE 
are also given. 
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a soft jellium density 



Pio 



(6) 



1 + exp((r - R(0))/a) 

where quadrupole and hexadecapole axial deformations are introduced through the jellium 
radius as 

R{6) =R [1 + 5 2 Y 20 (6) + 5 4 Y 4O (0)] . (7) 

The optimal deformation parameters 5 2 and £4 are determined by minimizing the total 
energy. 

We consider the clusters Na + X , Na + 5 , Na + 9 , Na + 7 , Na + 5 , Na + 5 , and Na + 19 , which, according 
to jellium estimates [EDI E2], I 



exhibit axial deformations. These clusters 



represent a wide size region and, as shown in Table I, cover prolate (Na + X , Na + 5 , Na + 7 , 
Na + 5 ) and oblate (Na + 9 , Na + 5 , Na + 19 ) shapes in ground states. Moreover, few of them (Na + 5 
and Na + 19 ) have shape isomers with a tiny energy deficit AE ~ 0.02 eV |2(J 22j and with 
quadrupole deformation of opposite sign with respect to the ground state. The largest 
sample Na + 19 has also a hexadecapole isomer. 
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Table I also shows the multipole moments 



4tt / drp {r)r x Y X0 ~ 5 j drp {r)r 2 
3 N e fr ' V 3 / dT P°W 

where A = 2,4 and Po( r ) is the ground state density of valence electrons. The dimensionless 
multipole moments f3\ are less model dependent because they are computed from expec- 
tation values. Thus they can serve for robust characterization of the deformation and for 
comparison between different models. The quantities 5\ in the jellium model @ coincide 
with the f3\ for small deformation. 

The optical response is calculated in the framework of the random phase approxima- 
tion (RPA). Full RPA for deformed systems is extremely involved. We employ a separable 
approximation of RPA where the residual interaction is expanded into a sum of separable 
terms |{Z0 |. The expansion employs local one-body operators Q\i p (r). Their structure is 



constructed as to map the response mean-field in RPA and the generating operators (to 
which response is explored) are chosen to cover the leading multipole operators. The expan- 
sion coefficients are computed self-consistently. It was shown that this procedure provides a 
sufficiently precise reproduction of the exact residual interaction [pOfl . I n particular, for the 
description of the SM, we use the basis of generating operators 

/ 2 ip = r 2+p (Y 21 (9) + Yl(e)), p = 0,2,4, 

Ui P = r 4+ P(Y 41 (9)+Yl(9)), p = 0,2. (9) 
The same set of operators was used for the description of Xp = 21 branch of the quadrupole 



plasmon p0| . The close similarity between scissors and quadrupole fields was discussed in 
The p = component of the input field /210 has a form of the external quadrupole field 
in the long-wave approximation. It generates the main piece of the separable interaction, 
peaked at the surface of the system. The next two quadrupole fields (p — 2 and 4) result in 
the separable operators Q2i P ( r ) peaked in more interior of the cluster. We include also hex- 
adecapole fields f^ p in order to account for coupling between quadrupole and hexadecapole 
modes, arising with the onset of deformation, especially in systems having both quadrupole 
and hexadecapole deformations. The set of the fields @ ensures good convergence of the 
separable expansion to exact results. Explicit expressions for Q\i p (r) are given in f2"0fl . 

We study the SM response in terms of photo-absorption. In axially deformed systems, 
the photo-absorption cross-section from the ground state to the excited state j = A 71 " of 



excitation energy uj is 



aiXX^gs^j) = J^ X+ } (^-f^\ < j\Ol\gs>\ 2 (10) 
v p>y A[(2A + 1)!!] 2 W 1 Jl v ' y 1 v ; 

where < jlO^lgs > is the reduced transition matrix element and is the operator of 
electric (X = E) or magnetic (X = M) transition. For the scissors mode, we have Of( = l x 
(in Bohr magnetons /ie). The selections rules are \i = A and (— 1) A = n for (X = E) or 
(— 1) A+1 = 7r for (X = M). The electric photo-absorption strength (10) will be used in 



Section |VII| for estimation of the competition between the scissors and low-energy electric 
excitations. 

Useful measures for the SM are provided by the sum rules 

5 m (Ml) = ^u;fS(Ml) i (11) 
j 

for m = —1,0, and 1. The ratios u = \J Si/ S_i or u = Si/ S provide a rough energy 
centroid of the low-energy mode. The 5_i is proportional to the paramagnetic susceptibility 
and Si to the integral photo-absorption cross section. The detailed values for the S m are 
obtained from explicit RPA results. A simple estimate for Si can be obtained using Eqs. 
(0) and (|): 

20 7 

Si(Ml) = ^Tc^Ml), = —NJlfil (12) 

3 Ts 

Sum rules do exist also for electric excitations where they can be reduced to the simple 



expressions 



Si(EX) = J2^\<j\ erXY ^\9s > | 2 (13) 

3 

h 2 e 2 



-A(2A + l) 2 X e < r lx - 1 > . 



IV. SCISSORS RESPONSE 

The Ml optical responses in light clusters are shown in Fig. 2. The present Kohn-Sham 
calculations yield results close to the ones obtained by using a deformed Woods Saxon 
potential Q. Only the excitation energies in prolate clusters are ~ 0.2 eV higher, because 
of the larger quadrupole deformations obtained in the self-consistent approach. The main 
characteristic of the Ml response is the occurrence of one or two prominent peaks below 1 
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TABLE II: Sum rules S m (in eV m ^1) calculated in the energy region 0-6 eV. The fractions for the 
region 0-1 eV are given in parenthesis. 
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eV. They are only slightly shifted from their unperturbed lph (particle-hole) spectrum. In 
the low-energy region, the lph [A^ = 1 + > spectrum is very dilute and, therefore, does not 
meet the conditions for Landau fragmentation or for pronounced coherent superpositions. 
Because of these features, the low-energy scissors strength can be associated to well defined 
lph transitions. In the Nilsson-Clemenger notation [Afn z A], they are [110] — > [101] in 
Na+, [211] -> [202] and [211] -> [200] in Na+ , [211] -> [220] in Na+, [321] -> [310] and 
[321] -»■ [312] in Na^ 7 , [310] -> [321] and [312] -»■ [321] in Na^ 5 . 

As shown in Fig. 2 (compare full and dashed lines), the residual interaction induces a 
rather moderate blue-shift that is much smaller than e.g. the shift for the electric dipole 
plasmon. The underlying physics is explained in Fig. 3 where the exchange-correlations and 
Coulomb contributions to the leading part of the separable operator Q 2 n ( r ) are presented 
for the case of Na^ 9 (the repulsive Coulomb is negative in this representation). It is seen 
that both contributions compensate each other to a large extent. The final outcome is a 
slight repulsive interaction responsible for the blueshift. It is worth noting that the balance 
between these two contributions is rather fragile and can be affected by different factors, e.g. 
triaxiality and detailed ionic structure, making even the sign of the net interaction uncertain 

Table II collects summed scissors strengths in wide- (0-6 eV) and low-energy (0-1 eV) 
regions, calculated within the SRPA. It is seen that the AJV = low-energy scissors mode, 
being mainly concentrated below 1 eV, contributes strongly to S-i and Sq. This justifies 
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using S-i and So for a rough estimation of the energy of the mode as u = a/ Si/S-i (or 
uj = S\/ Sq). Besides, since S-i is proportional to the paramagnetic susceptibility, this 
means that just the low-energy SM determines the van Vleck paramagnetism (see discussion 
in Section |VT| ). The high-energy part of the scissors strength (associated with AJ\f = 2 
transitions) contributes appreciably to the Si sum rule, i.e. to the total photo-absorption 
cross section. Table II also shows that the increase of the S m with cluster size is not 
monotonous. The fluctuations are caused by the predominantly lph character of the low- 
energy SM and by the difference in cluster deformations. NaJ 7 demonstrates especially 
strong Ml strength, while in oblate clusters Na^~ 9 and Na^ 5 the strength is relatively weak. 

Fig. 4 illustrates the coupling of the high-energy scissors branch to the quadrupole 
plasmon in Na^ 5 . The correlation between Ml and E2 peaks at ~ 3 eV is clearly seen. The 
coupling of electric and magnetic modes with the same quantum numbers A 71 " is a general 
feature of deformed finite quantum systems. It is well known, for example, in atomic nuclei 
(see, e.g. [|, |7|). 

The coupling of the SM with dipole and spin-dipole oscillations in clusters was discussed 
in detail in Ref. 0. In particular, it was shown that breaking the symmetry by the ionic 
lattice results in a weak coupling between these oscillations of opposite space parity. This 
feature may change for clusters deposited on a surface since symmetry breaking is much 
stronger there. 



V. EFFECTS OF SHAPE ISOMERS AND IONIC STRUCTURE 



The calculations |2(J, £2] show that light axially deformed clusters (N e < 40) have at 
the deformation-energy surface one distinct minimum corresponding to the ground-state 
deformation. The heavier clusters (40 < N e < 100 ) have usually two minima with opposite 
quadrupole deformations (prolate and oblate) and very close energies. The energy difference 
between the ground and the shape isomeric states is often less than 0.02 eV = 200 K. The 
number of shape isomers with a tiny energy deficit increases with cluster size. These isomers 
may have a variety of very different shapes. And significant amounts of isomers can be found 
in a thermal ensemblse and thus contribute to the SM, e.g., at room temperature. 

Fig. 5 compares the scissors modes built on the ground and isomeric states in Nas5 + . 
Both of them show a rich spectrum of low lying Ml states. But the spectra and typical 
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strengths of the two states differ substantially. Indeed, the Ml spectra are generated by 
several lph Ml transitions between the levels lying in the vicinity of the Fermi surface. 
Prolate and oblate shapes yield different sequences of the single-particle levels just at the 
Fermi surface. As a result, some of the lph Ml transitions which are significant in the prolate 
case are transformed to lpp or lhh transitions in the oblate case, thus strongly decreasing 
the strength. 

Fig. 5 compares also the SRPA results obtained with soft jellium model © and detailed 
ionic structure ||. In the case of the ionic structure, the global cluster deformations slightly 
deviate from those in the jellium case. It is seen that calculations with the ionic structure 
give somewhat different positions of the SM peaks with respect to the jellium case. Such 
a redistribution of the peaks is due to different energies of the single-particle levels in the 
ionic calculations. The photo-absorption response in the ionic case is stronger, which is also 
explained by the redistribution of the levels in the mean field. The level sequence in this 
case favors a few more lph transitions of Ml type. 

Altogether, the calculations demonstrate some important points concerning the SM in 
medium clusters: 

i) The mode is distributed over several lph transitions. 

ii) The single particle spectrum near the Fermi level is rather dense and so even small changes 
in the calculation scheme can redistribute visibly the spectrum and open (or close) some 
relevant lph transitions. 

iii) As was mentioned above, the ground and first isomeric states in medium clusters can 
be very close in energy. So, the scissors mode in free medium clusters should be considered 
as a statistical mixture of contributions from different shapes (predominantly of the ground 
state and first isomer). Such kind of the analysis (involving the deformation splitting, 
Landau fragmentation and contributions of shape isomers) was recently applied to explain 
experimental El optical response in deformed sodium clusters with 50 < N e < 60 P^| . 

Less ambiguous deformations will be provided by the SM in metal clusters deposited on 
insulating substrates |28|]. These clusters are oblate and their mean size and magnitude of 
the deformation can be well controlled. Such clusters seem to be the most promising systems 
for experimental search of the SM. 
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VI. MAGNETIC ANISOTROPY 



RPA calculations show that the SM energies and B(M1) strengths scale with the defor- 
mation <5 2 and the electron number N e basically according to the trends (Q) and (|j). Strong 
fluctuations, however, take place in small clusters M. They reflect the lph nature of the 
transitions and may affect the magnetic susceptibility. 

The total orbital magnetic susceptibility in clusters is the sum of Langevin diamagnetic 
and van Vleck paramagnetic terms [T], 



v _ ^,dia i para 

Xk-Xk +Xk > 



(14) 



where 



X k 



dia 



k ■ 



(15) 
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Xk 



< j\L k \Q > 
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the moment of inertia and by 



< j\L k \0 > 



j 



CO j 
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R 



N e <p 2 k > 



(16) 



(17) 



(18) 



its rigid-body value. The sum in (|T6| ) and ([[7]) runs over excited states. Further, k = x, y, z is 
the coordinate index. In axial systems, one has p 2 z = 2 < x 2 > and p 2 x y =< x 2 > + < z 2 >. 

Note that for k = x,y the operator entering in the matrix element in (|16|) is exactly the 
scissors generator. This makes evident that just the low-energy SM mainly contributes to 
X V xy a - Indeed Table II shows that contribution of the low-energy scissors mode to the value 
S-i ~ XT™ achieves 85 — 100%. So, just the SM determines the van Vleck paramagnetism. 

In the schematic model JI|, the moment of inertia comes out as the rigid-body value, so 
that 9 Xjy = 6 xy and, therefore, Xx% a = ~X d x % %i i- e - a complete compensation of dia- and 
para-magnetic terms in \x,y takes place. Due to axial symmetry, one also has x p z ara = 0- 
The total susceptibility becomes, therefore, strictly anisotropic fl|] 



X:,. 



Xy = 0, 



Xz 



dia 

Xz > 



(19) 



going from zero to diamagnetic values. 
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On the other hand, strong shell effects in some peculiar light clusters may alter apprecia- 
bly the above result. This is illustrated for Na^ 7 in Fig. 6. Because of very low excitation 
energy of the SM in this cluster (see an exceptionally large value of S_i for Na^~ 7 in Fig. 2 
and, also, in Table II, the paramagnetic susceptibility is enhanced considerably and is no 
longer balanced by the diamagnetic term. So, Na^~ 7 should be paramagnetic in x,y-directions 
and diamagnetic in z-direction. The cluster Na^ also hints this property. The magnetic 
moments in Na^ 7 are sufficiently large to allow for a measurement of the dia-para anisotropy. 
The observation of this effect would provide a strong (though indirect) evidence of the SM 
in clusters. 



VII. EXPERIMENTAL PERSPECTIVES 
A. General analysis 

As was mentioned in the Introduction, the SM is not yet observed experimentally in 
metal clusters. The search of the SM is hindered by several factors: 

1) The mode has very low photo- absorption cross section. Following our results, in sodium 
clusters with N e ~ 10 - 10 2 , a(Ml)/N e ~ 1(T 5 - 1(T 7 A 2 as compared to a(El)/N e ~ 2A 2 
for the dipole plasmon. Such a weak Ml signal is at the edge of the sensitivity of modern 
detectors. 

2) Eqs. (|T0|)-(|T3|) allow to estimate the ratios for maximal optical responses as 

a(El)/a(Ml) = 0.35 • 10 5 A~ 2 (^) 2 (20) 

$2 



and (for u E2 ~ 3 eV §, §) 

2 

a{E2)/a{Ml) = 0.82 ■ 10~ 2 A- 4 (^) 2 iV 2 / 3 . (21) 

This gives for deformed (82 = 0.2) sodium (r s = 2.lA) clusters a(El)/a(Ml) ~ 4 • 10 6 and 
a(E2)/a(Ml) ~ 10 2 (the latter for N e = 125). So, the SM suffers from the competition 
with El and E2 strengths. Besides, the competition with E2 increases with cluster size. The 
above estimates compare the SM optical response with the maximal responses of El and 
E2 plasmons. As is shown in the next subsection, the competition is much weaker in the 
low-energy region where the SM has its stronghold, but one has still to be aware of large 
amounts of El strength. 
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3) The low-energy scissors mode lies in the infrared region where commonly used detectors 
for a visible light are not efficient enough. 

4) The SM energy decreases with cluster size and reaches the region of phonon excitations, 
~ O.leV, in very big clusters. 

In spite of these hindrances, clusters offer enough opportunities to look for the optimal 
conditions for observing the SM: one may change cluster size and (or) deformation, use 
different metals, choose between free and supported clusters, etc.. The macroscopic estimates 
(!)"(!) an d the present jellium RPA calculations may serve as a guide. Such analysis is briefly 
given below. To this purpose, we consider, as a typical case, sodium clusters with a moderate 
deformation 62 = 0.2. 

Larger clusters help because the Ml photo-absorption cross section grows linearly with 
the size. Already in deformed clusters with N e = 10 4 — 10 6 the scissors signal should be 
detectable. Unfortunately, the corresponding energy uo ~ 0.1 — 0.01 eV approaches or even 
covers the region of phonon excitations. Thus the best compromise for free sodium clusters is 
achieved at N e ~ 10 3 . Besides that, these clusters are large enough to ensure the dominance 
of the orbital scissors mode over spin Ml excitations. 

We may also look for clusters of larger density and deformation to increase the energy 
and the strength of the mode. For example, Li clusters (r s = 1.7 A) allow to increase both 
the energy and optical response by a factor of ~ 1.5 with respect to Na (r s = 2.1 A). Further 
improvement may be achieved with Ag (r s = 1.6 A), Mg (r s = 1.4 A), or Al (r s = 1.1 A). 
In any case, highly deformed clusters are welcome because a (Ml) ~ b\ and u> ~ 62- 

In general, free clusters seem not suitable for observing the SM. They are well mass- 
separated only up to sizes of hundred atoms. However, these systems have a weak Ml 
signal. As for heavier clusters, they suffer from a poor mass separation and are expected 
to be weakly deformed [26]. Moreover, as was discussed in Section [V], their Ml signal is a 
statistical mix of contributions of different shapes given by the ground state and isomers. 

Deposited clusters look more promising. One can adopt techniques that allow to get 
oblate clusters (Na and Ag) on dielectric surfaces |28j and, more remarkably, to monitor 
their size and deformation. Clusters with 10 2 — 10 6 atoms can be used for this aim. In this 
way, one can obtain supported clusters of a desired size and shape. High density of clusters 
on the surface gives good statistics in the measurements. Last but not least, monitoring 
the shape of clusters gives the chance to use the trends a (Ml) ~ b\ and lumi ~ $2 to 
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distinguish the SM from the El and E2 signals. To this end, one has to irradiate clusters of 
the same size but of different deformations. Then the low-energy El and E2 cross sections 
are also affected but, hopefully, in an irregular fashion (thus averaging out) or, at least, in 
a way different from the trends mentioned above. The cross sections measured at different 
deformations can be mutually subtracted to extract a small useful signal. For the analysis of 
such experiments, we need accurate theoretical estimates accounting for all the main factors 
and, in particular, the influence of the surface on the electric and scissors plasmons. Such 
work is in progress. 

The competition with El and E2 modes may be bypassed partly by resorting to specific 
reactions and techniques capable of hindering if not suppressing the El and E2 channels. 
We will discuss this issue for three relevant reactions: photo-absorption, resonance Raman 
scattering, and inelastic electron scattering. 



B. Photo-absorption 

The SM lies much lower than the dipole and quadrupole plasmon modes (0.2-1.0 eV 
against 2.5-3.0 eV for the dipole and 2-4 eV for the quadrupole fl9], |2(|). In principle, one 



could exploit this energy separation to focus on the SM. However, its Ml photo-absorption 
cross section is still very small as compared to the El strength and can be masked by the 
tail of close El modes. This is illustrated in Fig. 7 which compares El, E2 and Ml photo- 
absorption strengths calculated within SRPA approach pDJ in axially deformed Na+ 19 (ground 



state). The El strength, though extremely small compared to the dipole plasmon peak, is 
still strong enough to mask the SM. The competition with the quadrupole E2 strength in 
Na^g is negligible (though it can be stronger in heavier clusters). 

The competition with the El channel can be decreased by using deformed (oblate) clusters 
with N ~ 10 3 — 10 4 atoms, supported on dielectric surfaces, and adopting infrared techniques 
with polarized light. For the SM, the magnetic field of the incoming light should be parallel 
to the surface. This can be done with both s- and p-polarizations (with the electrical field 
to be perpendicular and parallel to the incident plane, respectively). The s-polarized light 
under an angle ~ 0° with the surface normal suppresses the \l — El branch without 
weakening the SM. The p-polarized light at angle ~ 90° suppresses the fi = 1 El branch, 
again without weakening the SM. The second variant seems to be preferable since the n = 1 
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branch is generally stronger (see, e.g., Fig. 7. 



C. Resonant Raman scattering 

In resonance Raman scattering (inelastic scattering of polarized photons), the SM may be 
populated through the de-excitation of the dipole plasmon excited by the incoming photon. 
The difference between energies of incoming and outgoing photons gives the energy of the 
mode. One must look at El decays, which represent the leading channel and may populate 
in deformed systems Ml as well as E2 states. With respect to photo-absorption, resonance 
Raman scattering has the advantage of excluding the competition with the strong El mode 
and to work in the visible light region. An experiment of this kind was carried out in 
supported Ag clusters |29j . The analysis of the data suggested, however, that the E2 rather 
than the SM was populated through the El decay. On the other hand, it is not simple to 
separate clearly the E2 from the Ml modes which are mixed in deformed systems. 



D. Inelastic electron scattering 

Inelastic electron scattering is also a promising technique. At backward scattering angles, 
the Ml signal should prevail over the electric ones. This property enabled the discovery of the 
scissors mode in deformed atomic nuclei 0. Polarization of the electrons provides additional 
possibilities to extract the magnetic response. In nuclear physics, electrons with very high 
energies (tens and hundreds MeV) are used, which favors creation of well collimated, intense 
and monochromatic electron beams. Moreover, at so high energies, it is easier to enhance 
the relative contribution of the transverse form factor. 

The low-energy, intense and monochromatic electron beams, suitable for our purpose, as 
well as the angular resolved techniques are now available |J{J. To our knowledge, however, 
inelastic electron scattering was used so far only to observe the dipole plasmon. Theoretical 
estimates |3l| show that at certain scattering angles the dipole contribution can be well 



suppressed in favor of the quadrupole plasmon, which, therefore, can be observed. From 
the latter analysis, we may infer that also the SM may be extracted. Indeed, in the low- 
energy region the SM can be stronger excited in the photo-absorption than the E2 mode. 
The analogy with atomic nuclei suggests that it should be favored even more in electron 
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scattering at backward angles. 



VIII. CONCLUSIONS 

The scissors mode (SM) in axially deformed sodium clusters has been studied within a 
self-consistent RPA approach based on the soft jellium model for the ionic background and 
density-functional theory for the electrons. The calculations show that the distribution of 
the low-energy scissors strength is dominated by lph states with only small shifts by the 
residual interaction. The SM thus reflects almost directly the density of lph Ml transitions 
near the Fermi surface and exhibits an extreme sensitivity to the actual deformation. The 
SM is rather unambiguous in light clusters (N e < 40) while it should be a statistical mix of 
contributions from different shapes in heavier clusters. Furthermore, the high-energy branch 
of the SM in deformed systems is coupled with the quadrupole plasmon. 

The calculations show that the SM determines the Van Vleck paramagnetism and results 
in a strong anisotropy of the magnetic susceptibility. Moreover, due to quantum shell effects, 
some clusters (Na^ 7 ) demonstrate dia-para anisotropy when the axial cluster is diamagnetic 
along the symmetry axis and paramagnetic along other two axes. The magnetic moments 
are large enough to be measured. The experimental observation of the dia-para anisotropy, 
being interesting itself, could also serve as an indirect evidence of the SM in clusters. 

A general analysis of possible routes for experimental observation of the scissors mode was 
done. The photo-absorption with a polarized light, inelastic electron scattering and Raman 
scattering were considered. The former two reactions seem to be more promising. Moreover, 
deposited clusters favor the detection of the SM. Its features need yet to be worked out in 
detail. Work along that direction is in progress. 
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FIGURE CAPTIONS 



Figure 1: Macroscopic view of scissors mode : rigid rotation (left), and rotation 
within a rigid surface [I| (right). 

Figure 2: Photo-absorption cross-section for the scissors mode, weighted by the Lorentz 
function with the averaging parameter 0.25 eV. The responses with (solid curve) and without 
(dashed curve) the residual interaction are presented. 

Figure 3: Radial profile of the residual interaction in Na^ 9 . Exchange-correlation (long- 
dashed curve) and Coulomb (shot-dashed curve) contributions are given together with their 
sum (solid curve). 

Figure 4: Coupling between the high-energy scissors branch and E2 quadrupole plasmon 
in prolate Na^ 5 . Ml and E2 photo-absorptions are weighted by the Lorentz function with 
the averaging parameter 0.25 eV. 

Figure 5: The strength distribution of the scissors modes in Na^" 5 , built on the prolate 
ground state (upper panel) and oblate first isomer (lower panel). Jellium (solid curve) and 
ionic (dashed curve) SRPA results are compared. 

Figure 6: Normalized diamagnetic, paramagnetic and summed moments /i = \B {B = 
AT) in axial deformed clusters Na^, Na^~ 5 , Na+ 9 , Na^~ 7 , Na^~ 5 , Na^ 5 and Na^ 19 . 

Figure 7: Photo-absorption El, E2 and scissors cross sections in Na^g weighted by the 
Lorentz function with the averaging parameter 0.25 eV. The El \l — and 1 branches are 
given separately to show their relative contributions to the competition with the scissors 
mode. 
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